/*******************************************************************************
																				
	DESCRIPTION:  	This file produces panel B of Figure 2.
	
*******************************************************************************/

clear all
global id_code 112
	
*******************************************************************************
 * Produce histograms comparing the 2006 and 2009 model
*******************************************************************************

* Combine the predictions	
use "${data}/003_MainWithEnsemblePred_Full_2006.dta", clear
rename p_emplAft6M_0M_In p_emplAft6M_0M_In_2006

append using "${data}/003_MainWithEnsemblePred_Full_2009.dta"
rename p_emplAft6M_0M_In p_emplAft6M_0M_In_2009

* Make the graph:
twoway ///
	(kdensity p_emplAft6M_0M_In_2006 , recast(area)  fcolor(ebblue%30) lcolor(ebblue)) ///
	(kdensity p_emplAft6M_0M_In_2009 , fcolor(orange_red%30) lcolor(orange_red) recast(area)), ///
	graphregion(color(white)) plotregion(margin(b=0 l=0))						///
	xtitle("6-Month Job-Finding Probability at Start of Spell (Baseline Model)") ytitle("Density") 			///
	xscale(titlegap(2)) yscale(titlegap(2)) ///
	xlabel(0(0.1)1) ylabel(0(1)3, angle(0) format(%9.0g)) ///
	legend(cols(2) order(1 "2006" 2 "2009") region(lwidth(none)))

graph export "${output}/${id_code}_Predicted6MJFR_Distribution_0M_Full_2006vs2009_kdensity.pdf", as(pdf) replace
